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^ \ Abstract: Cosmic ray events above 10^° eV are on the verge of confronting 

fundamental particle physics. The neutrino is the only candidate primary 

■ among established particles capable of crossing 100 Mpc intergalactic dis- 

2 ■ tances unimpeded. The magnitude of vN cross sections indicated by events, 

O ■ plus consistency with the Standard Model at low-energy, point to new physics 

O 
O 
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of massive spin-2 exchange. In models based on extra dimensions, we find 
that the vN cross section rises to typical hadronic values of between 1 and 



p i' 100 mb at energies above 10 eV. Our calculations take into account con 



straints of unitarity. We conclude that air-showers observed with energies 



D . above 10^^ eV are consistent with neutrino primaries and extra-dimension 

^ ' models. An upper bound of 1-10 TeV on the mass scale at which graviton 



exchange becomes strong in current Kaluza-Klein models follows. 

1 Introduction 



The energy of extra-galactic proton cosmic rays should not exceed the GZK 
bound [^. The bound, about 10^^ eV, is based on the known interactions of 
nucleon primaries with the photon background of intergalactic space. The 
GZK bound is tantamount to an upper limit on cosmic ray energies, inas- 
much as nuclei and photons have lower energy cutoffs . Yet an experimental 
puzzle exists, as evidence for air shower events with energies above the GZK 
bound has steadily accumulated over the last 35 years There seem to 
be inadequate sources nearby to account for such events, and the sources 
are almost certainly extragalactic 0]. A number of showers with energies 
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reliably determined to be above 10^° eV have been observed in recent years 
[0], deepening the puzzle. 

A completely satisfactory explanation of the so-called GZK-violating 
events is still lacking. Models have been constructed to explain the puzzle by 
invoking "conventional" extragalactic sources of proton UHE acceleration, 
as reviewed recently in Ref.0. Other models introduce exotic primaries such 
as magnetic monopoles or exotic sources such as unstable superheavy relic 
particles 0, or appeal to topological defects 0. Except for monopoles all of 
the proposed sources must be within 50-100 Mpc to evade the GKZ bound, 
a requirement which is difficult to satisfy. 

The 10^° eV events potentially pose a confrontation between observa- 
tion and fundamental particle physics. Except for the neutrino, there are 
no candidates among established elementary particles that could cross the 
requisite intergalactic distances of about 100 Mpc or more [110, ITT| . The neu- 



trino would be a candidate for the events, if its interaction cross section were 
large enough, but this requires physics beyond the Standard Model. The 
neutrino-nucleon total cross section atot is the crucial issue: Flux estimates 
of UHE neutrinos produced by extra-galactic sources and GZK-attenuated 
nucleons and nuclei vary widely, but suffice to account for the shower rates 
observed. Very significantly, there is a hint of correlations in direction of the 



events, both with one another and candidate sources |T^. Event-pointing 
toward distant sources, if confirmed, would require a neutral, long-lived pri- 
mary, reducing the possibilities for practical purposes to the neutrino plus 
new physics to explain the cross section. 

Current understanding of the UHE Standard Model atot is based on 
small-x QCD evolution and W^,Z exchange physics [0. This physics is 
extremely well understood and has been directly tested up to s = 10^ GeV^ 
with recent HERA-based updates [0]. These calculations are then extrap- 
olated to the region of 10^° eV primary energy, leading to cross sections in 
the range 10~^ — 10~^ mb, far too small to explain the GZK-violating air 
shower events. The observationally indicated cross section is completely out 
of reach of the extrapolations of W- and Z- exchange mechanisms. 

Since the neutrino-nucleon cross section at 10^° eV has never been directly 
measured, it is quite reasonable to surmise that new physical processes may 
be at work. Total cross sections at high energies are dominated by charac- 
teristics of the t-channel exchanges. The growth of atot with energy, in turn, 
is directly correlated with the spin of exchanged particles. Exchange of new 



2 



(W- or Z-like) massive vector bosons would produce atot growing at the same 
rate as the standard one, faihng to explain the puzzle. If the data indicates 
a more rapid growth with energy, one is forced to consider higher spin, with 
the next logical possibility being massive spin-2 exchange. We reiterate that 
this deduction is data-driven; if data indicates (a) corrrelations with source 
directions, and (b) atot in the mb and above range, there are few options 
other than neutrinos interacting by massive spin-2 exchanges. 

Recent theoretical progress has opened up the fascinating possibility of 
massive spin-2, t-channel exchange in the context of large "extra" dimensions 



rH|, while the fundamental scale can be related to a string scale of order 



several TeV |]T6|. In this context the Kaluza-Klein (KK) excitations of the 
graviton act like a tower of massive spin-2 particle exchanges . We will show 
that large UHE neutrino cross sections, sufficient to generate the observed 
showers, are a generic feature of this developing framework [T^ . At the same 
time the new contributions to atot at energies below center of mass energy 
a/s of 500 GeV is several orders of magnitude below the Standard Model 
component. In fact, the new physics we propose to explain the puzzle of the 
G Z K-violaiing events is consistent with all known experimental limits. 



2 The uN Cross Section with Massive Spin-2 
Exchange 

The low-energy, 4-Fermi interaction total cross section, atot, grows like over 
many decades of energy. Perturbative unitarity implies that at an invariant 
cm energy ^/s large compared to the exchange mass mw, the growth rate 
slows to at most a logarithmic energy dependence. The shift from power- 
law to logarithmic growth is seen to occur in the Standard Model. There 
is a second effect, that above 100 TeV the total number of targets (quark- 
antiquark pairs) grows roughly like {E^)^''^. This fractional power, in turn, 
leads to a formula for atot = 1-2 x 10~^mfe(£'j,/10^^eV)*^'^ as a reasonable 
approximation to the Standard Model calculation [|1^ . 



Exchange of additional spin-1 bosons cannot produce faster growth with 
energy than just described. However, a massive spin-2 exchange grows quite 
quickly with energy on very general grounds. A dimensionless spin-2 field 
gets its couplings from derivatives, which translate to factors of energy. Thus 
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the naive cross section grows like E'^, in the "low-energy" regime. 

These general features are exemplified in the Feynman rules for this 
regime developed by several groups |TB[. To be consistent with the liter- 



ature, we will describe the interaction as "graviton" exchange, implying the 
standard picture of a tower of spin-2 KK modes. The parton level v gluon 
differential cross section is given by. 



da^3 ttA^ u 



di 2M% P 



2u^ + AuH + 3mP + 



(1) 



Here Ms is the cutoff on the graviton mass, and A is the effective coupling at 
the scale Ms that cuts off the graviton KK mode summation. The magnitude 
of parameter A has been lumped into the scale parameter M5, hence A = ±1 
for our purposes ||19|. In Eqs. (1) and (2) we take —t <C M|, which leads 
to the simple factor 1/M|. This suffices for our extrapolation, but the full 
t dependence is used in the partial wave projections to check the unitarity 
constraint (discussed momentarily). The corresponding parton level v quark 
differential cross section is given by. 



— ^ = — __ 32m^ + UuH + A2uh^ + im^ + 
dt 32Mg s2 L 



(2) 



We include the contribution of the two valence quarks as well as the m, d and 
s sea quarks. The Z-graviton interference terms are included with negative 
A, though their contribution is very small compared to other terms. The 
negative sign gives a slight enhancement for the final result of a tot- Collider 
physics and astrophysics constrain the effective scale M5 to be above 1 TeV, 
with lower number of dimensions leading to stronger constraints. 



2.1 Unitarity 

The complete theory of massive KK modes is not yet developed, making 
it impossible to know the exact cross sections at asymptotic energies. The 
situation is analogous to the case of the 4-Fermi theory before the Standard 
Model. By observing the growth of atot it was possible to deduce a massive 
vector exchange long before a consistent theory existed. In much the same 
way, present data indicate a spin-2 exchange while the analogous complete 
"standard model" of gravitons does not yet exist. Unlike the electroweak 
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case, in either the data- driven or extra-dimensions scenario we must face a 
strongly interacting, non-perturbative problem in the high energy, s ^ M|. 
Perturbative unitarity breaks down as a host of new channels opens up in 
that regime. The low-energy effective theory remains an accurate description 
within a particular domain of consistency. Extrapolation of the 4-Fermi 
predictions to higher energies is possible by matching the consistent, low- 
energy description with the asymptotic demands of unitarity. Similarly, we 
resolve the difficulties of massive graviton exchange in the high energy regime 
by matching the y/s < Ms predictions, where the perturbative calculation is 
under control, to the '^/s ^ Ms non-perturbative regime. 

We proceed by first evaluating the theory's partial wave amplitudes to 
find the highest energy where the low energy effective theory is applicable. 
Taking the case u + q ^ u + q, and including the full dependence of the 
propagator, we find that the unitarity bound on the J = projection of the 



helicity amplitude T++ |^ gives the strongest bound. For example, with 
the number of extra dimensions n = 2 we find s/ s < 1.1 Ms-, while with n = 4 
we find \/s < 2.OM5. As mentioned earlier, the most attractive value of Ms 
is in the TeV range. 

The invariant energies of the highest energy cosmic rays are approximately 
1000 units of the scale Ms ~ 1 TeV, well beyond the low-energy regime. A 
phenomenological prescription consistent with unitarity is clearly necessary 
to extrapolate the low energy amplitudes. We now turn to describing and 
motivating three different asymptotic forms that span reasonable possibili- 
ties: log{s), and growth of <7tot{s). There is no guarantee a priori that 
any should extrapolate from low to high energy with Ms ~ 1 TeV and pro- 
duce hadronic-size cross sections at 10^° eV. As we shall see, surprisingly, 
they all do! 

As a first version of an extrapolation model, we use a well known result 
from general features of local quantum field theory. The Froissart bound 
[PHI, reflecting the unitarity constraint on cross sections from exchange of 
massive particles, dictates that atot grows no more rapidly than (log(s/M|))^. 
The bound is an asymptotic one, and strictly speaking incapable of limiting 
behavior at any finite energy; moreover, the bound is probably violated in 
the case of graviton exchange. It is quite conservative to use the Froissart 
bound, with its mild logarithmic growth in s, as a first test case. In terms of 
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the differential cross section, we have at high energy 

da const , 
di ~* iMl 



\og{s/Ml) 



(3) 



We then propose the following interpolating formula which reproduces Eq. 
(1) in the low energy limit and Eq. (3) in the high energy limit 



da^3 
di 



u 



2Ml{Ml + sy{Ml-t)s^ 
'l + eiog(l + s7M^ ' 



(4) 



The i/-quark parton level cross section is similarly extrapolated to high en- 
ergies. We have introduced the parameter ^, which we will allow to vary be- 
tween 1 and 10. It cannot be much larger than 10 since then the low energy 
cross section gets modified, violating consistency. We use these parton-level 
expression to calculate a tot- We convolute the parton-level cross section with 
CTEQ 4.6 parton distribution functions, which give a continued growth in 
the UHE regime from the small-x effect. As pointed out previously, the cross 
section can be expected to grow to "strong interaction" magnitudes, where 
parton coalescence and string effects ultimately come into play. 

A different constraint for unitarizing would be growth. Regge theory 
would suggest growth for spin-2 exchange at small, fixed t, which (in 
fact) occurs in this theory when the t values are restricted self-consistently. 
Thus the use of growth makes a comparatively mild alteration of the 
perturbative predictions and follows from eikonal unitarization of Reggeized 
graviton exchange [^. Another unitarization procedure indicates a linear 



growth in s I^S], which represents a case intermediate to the other two. These 
cases serve to establish a fair range of possibilities for models of unitarization. 
For the and models, the extrapolation form of Eq. (1) we choose is 



da^s 
di 



u r 



2u^ + AuH + 3m? + F 



(5) 



where p = 1, for s^, s^. A similar extrapolation is applied to the cross 
section formula for the v quark parton case, Eq. (2). 

Note that we use the detailed perturbative low energy calculations that 
follow from |T8[ to anchor the low energy end. As we discuss below, we also 
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explore the allowed range of Ms values, corresponding to a range of numbers 
of extra large dimensions. As far as we know, these consistency features in 
the present context have not previously appeared in the literature. 

An essential feature of spin-2 exchange, complementary to the growth of 
atot with energy in the U HE region, is suppression of observable effects in the 
low energy regime of existing data. Turning the problem around to a "data 
driven" view, the energy dependence of the new physics must be so strong 
that the miUibarn-scale total cross sections at -y^s > 10^ TeV are suppressed 
well below the Standard Model values below 1 TeV. This also follows in our 
approach. 

3 Results 

Results of the calculations based on our models (Eqs. 4 and 5) are given 
in Fig. 1. For the ln{s) model, shown by the dotted line, we find that for 
incident neutrino energy of order 10^^ GeV, Otat is roughly 0.5 mb, with the 
effective cutoff scale A4s = 1 TeV and the choice ^ = 10. This cross section 
is remarkably consistent with the low end of the range required to explain 
the GZK violating events. Looking at higher values of Ms, we find that the 
Utot falls so steeply as a function of Ms that Ms ~ 1 TeV is required for this 
model to be a viable explanation of cosmic ray events with £■ > 10^^ eV. 

The results of the calculations with the model are shown in the dash- 
dotted line {j3 = 0.1) and short-dashed line {f3 = 1.0). The value of Ms 
is fixed at 1 TeV; (3 values of 0.1 and 1.0 are chosen to illustrate that this 
intermediate growth model easily produces hadronic size cross sections above 
the GZK cutoff. We find that raising the cutoff Ms much above 2 TeV 
suppresses the cross section well below the range required for the GZK events, 
independently of the value of f3. 

The long-dashed curve (Fig. 1) shows the case with growth for values 
Ms — 3 TeV and P — 1.0. The growth of atot to tens of millibarns at the 
highest energies shows the viabihty of this model even at Ms values above 1 
TeV. 

The key feature we emphasize about Fig. 1 is that atot somewhere in the 
range 1-100 mb is obtained at the highest energies for reasonable parameter 
choices for all three unitarization models. As a significant corollary, we find 
that if Ms is larger than 1-10 TeV (depending on the model), then the 
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low scale gravity models fail to increase atot by enough to explain the GZK- 
violating events. Thus the analysis serves to put an upper hound on the cutoff 
massQ. 

For completeness let us note that the interaction cross section with CMB 
photons is completely negligible because the cm-energy is about 10~^^ that 
of a proton target. On the other hand, the mean free path of UHE neutrinos 
with El, > 10^^ eV and m6-scale cross sections is of the order of a few kilo- 
meters in the upper atmosphere. This is still much larger than the mean free 
path of other particles such as protons at the same energy. If the ultrahigh 
energy showers are indeed initiated by neutrinos, the log{s) and models 
predict considerable vertical spread in the point of origin of the shower. The 
model is capable of much larger atot- 

Experimental extraction of cross sections is complicated by shower devel- 
opment fluctuations, but it would be useful to make all possible efforts to 
bound the mean free path of the primaries responsible for GZ/T-violating 
events. There is a further, interesting variant of the "double-bang" v-r signal 
2^ : secondary events may come from a neutrino dumping part of its energy 



in a first collision, undergoing a second collision a few kilometers later. This 
process may serve to separate primaries with m6-scale cross sections from 
protons or nuclei, which have much shorter mean free paths. 

3.1 Signatures of Massive spin-2 Exchange at Km-Scale 
Neutrino Detectors 

It is also interesting to consider alternate signatures of the v — N neutral- 
current cross section exceeding the Standard Model. The cosmic neutrino 
detectors, AMANDA [^] and RICE for example, are expected to explore 



the TeV-10 PeV energy regime. The cross sections we find exceeding those 
of the Standard Model might be tested in these experiments. Of particular 
interest is the angular distribution of events. The diffuse background cosmic 
ray neutrino flux is expected to be isotropic. We can, therefore, measure the 
deviations of atot from SM predictions by measuring the ratio of upward- 
to downward- moving events. This ratio is plotted in Fig. 2 as a function 
of the incident neutrino energy. The plots show a few representative choices 



^Our requirement is consistent with estimates of lower bounds as discussed in ^J], for 
example. 
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of Ms and extrapolation parameter f3, using the model for illustration. 
The up/ down ratio starts to deviate very strongly from the SM value for an 
incident neutrino energy greater than about 5 PeV. Beyond about 5 PeV the 
ratio falls very sharply to zero. This in principle can be measured at RICE 



[ ^ , which is sensitive to precisely the energy regime of 100 TeV to 100 
PeV. The event rate in this energy region is also expected to be significant. 
Note that only the neutral-current events are affected. A more detailed, 
but also attractive extension of this technique is UHE Earth tomography, 
recently shown practicable with existing flux estimates, and also capable of 



measuring atot indirectly The graviton-exchange predictions are not 

especially sensitive to the precise value of (3 in this region, but do depend 
strongly on the cutoff scale Ms- Large scale detectors such as RICE will be 
able to explore the range of Ms < 2 TeV. 

We conclude by reiterating that the highest energy cosmic rays are on 
the verge of confronting fundamental particle physics. Exciting projects un- 
derway, including AGASA, Hi- RES, and AUGER should be able to collect 
enough data to resolve the issues. The puzzle of GZK violating events can 
be experimentally resolved by establishing neutral primaries via angular cor- 
relations. That feature would imply new physics of neutrinos, with the cross 
section then indicating massive spin-2 exchange rather directly. Experimen- 
tally bounding the primary interaction cross section below that of protons is 
another viable strategy in the and log{s) scenarios. 

Current limits on KK modes of extra dimension models predict suffi- 
ciently large u — N cross sections to produce such interactions, and are 
consistent with observation of events in the upper atmosphere at primary 
energies greater than lO^'' — 10^^ eV. Depending on the observational out- 
come, then, the subject matter could become very important, and further 
development is well warranted. 
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Figure 1: The vN cross section in the Standard Model iSM) compared 
to a theory with large extra dimensions and three different models for the 
unitarity extrapolation between perturbative to non- perturbative regimes. 
The dotted line shows the log{s) growth case with Ms = 1 TeV and ^ = 
10. The short dashed and dash-dotted lines show growth with Ms = 1 
TeV and P — 1 and 0.1 respectively. The long dashed line shows growth 
with Ms — 3 TeV and P — 1. The contribution from massive graviton 
exchange is negligible at low energies but rises above the SM contribution 
when a/s > Ms, reaching typical hadronic cross sections at incident neutrino 
energies in the range 5 x 10^^ to 5 x 10^° GeV. The HERA data point is 
shown for comparison. The approximate minimum value required for uN 
cross-section, a — 1 mb, is indicated by the horizontal straight line. 
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Figure 2: The ratio of upward- to downward-moving events at UHE cosmic 
ray neutrino detectors as a function of the incident neutrino energy. The 
sohd curve corresponds to the SM prediction. Other curves include the 
contribution due to low scale gravity models with cross section growth 
for several different representative choices of the cutoff parameter Ms and 
the dimensionless extrapolation parameter j3 = 1. Since the extrapolation is 
small, the result is not especially sensitive to the precise value of (3. 
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